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Magnetic Reconnection in Different Environments: 
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15 . 1 . INTRODUCTION 
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Magnetic reconnection is broadly recognized as a key 
plasma transport and energy conversion process in space 
plasmas. Ranging from solar activity [e.g., Cargill and 
Klimckuk, 1991; Ant lochos ei al,, 1999] tliroxigh inter pi an* 
etary space [e,g,, Gosling et al, 1995, 2005] to the bound- 
aries [e,g,, Paschmann et al^ 1979; Sonnerup et aL, 1981] 
and interiors of planetary magnetospheres [e.g„ Hones, 
1977], magnetic reconnection facilitates highly time 
dependent, eruptive behavior and it transports mass, 
momentum, and energy into regions featuring otherwise 
nearly impenetrable boundary layera — 

A key feature of magnetic reconnection is that it 
primarily converts magnetic energy in the two inflow 
regions into internal and kinetic energy. The efTicacy of 
this energy conversion process depends strongly on 
plasma parameters in the two inflow regions [e,g,, Aunai 
et a/,2011]. For example, high plasma “ 2p^iEP [Sw/ 
et a/„ 2003] or low Alfven speed [Hesse et al, 1996] in the 
inflow slows down the reconnection process. 

In Earth’s magnetosphere, magnetic reconnection 
operates both at the magnetopause, the boundary bet- 
ween solar- wind*dominated plasma and the region dom- 
inated by Earth’s magnetic field, and on the nightside, 
within Earth’s magnetotail. In the former situation, 
reconnection merges magnetic fields and plasmas of solar 
origin on one side with magnetospheric magnetic field 
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and plasma on the other in a usually very asymmetric 
geometry [e,g,, Mozer et al 2008]. Furthermore, recon- 
necting magnetic fields are typically not coplanar, 
implying, in the reconnection plane> the existence of an 
out-of-plane magnetic field. On the nightside, however, 
reconnection is believed to involve considerably more 
symmetric environments, with usually relatively small 
guide magnetic fields. It is expected that similar situations 
prevail at the other magnetized planets. Given the varied 
nature of the environments magnetic reconnection oper- 
ates in, it is therefore of great interest to understand how 
reconnection adjusts to plasma parameters in both inflow 
regiona 

Owing to its importance, reconnection has therefore 
been studied by means of spacecraft observations [e,g,, 
Eastwood et al., 2007, 2010; Faschmann, 2008], by analytic 
theory [Cassak and S/iay, 2007; Birn et al., 2011; Hesse 
et al., 2011], and by numerical modeling. Recent work has 
indicated that reconnection changes substantially if the 
inflow sides are unequal et a/,, 2003; Mozer and 

Pritchett , 2011; Hesse et a/„ 2013] or if a guide magnetic 
field is being introduced [e.g., Hesse etal., 2011], Kinetic, 
hybrid, and HalfMHD simulations consistently exhibit a 
reconnection rate much larger than typical Sweet-Pa rker 
rates [Birn et ah, 200 1|; however, recent evidence suggests 
that island formation may impact reconnection rates in 
kinetic models [Daughton et aL, 2006, 2009, Karimabadi 
et a/,, 2007], and it may play a critical role in enabling 
fast reconnection even in MHD systems [Huang and 
Bhattacharjee, 2013]. 

This chapter therefore addresses two important aspects 
of magnetic reconnection: time-dependent rates affected 
by island formation and the preferred direction of the 
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iiaconiiectioii line if merging magnetic field components 
are not antiparallei The former presents new results, 
whereas the latter will summarize recent findings 
et al, 20 1 3]. The two aspects of the analysis will be prefaced 
by a basic discussion of the impact magnetic reconnection 
has on energy transport across topological boundaries. 

15.2. ENERGY TRANSPORT 

While it is widely known that magnetic reconnection 
facilitates mass, momentum, and energy transport in 
plasmas, it is illustrative to consider an analytical 
approach to describing the efficacy of reconnection as a 
transport mechanism. For simplicity, we consider a 
simple, two-dimensional geometry, where an X-shaped 
reconnection region forms the intersection of two sepa- 
ratrices. The geometry is depicted in Figure 15.1 . 

Here, V denotes the volume enclosed by the two sepa- 
ratrices, ^ V its boundary (the separatrices), n the local 
boundary normal, v the plasma velocity, and the 
velocity of the separatrix surface itself Note that the 
normal here is represented by the blue arrow in Figure 15.1 
and must not be confused with the normal direction to 
the reconnecting current sheet. 

Integrating the MHD energy equation over the volxime 
V then leads to 


V ^ p 
p— + 

2 7-1 
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— \dv 

dti 

= -T^x£ + — f^-vl + p^(ii,-v) 

iv L 2 ^ 

7-1 

Here we have assumed the generic relation 



two-dimensional system, we can represent the magnetic 
field through a flux function A such that 


B = VAxe,, 


(15.3) 


where is the unit vector in the direction of the ignor- 
able coordinate y. In the Coulomb gauge, the time evolu- 
tion of the magnetic field can then be described by 


— +v-Vv4 = -ii„. 

dt " 



(15.4) 


At the X point proper, the gradient of ^ vanishes. At 
the same time, the numerical value A of the flux function 

JC 

at the X point is constant along the evolving separatrix 
surfaces, implying 


dt 




^ lx point 


(15.5) 


Away from the diffusion region, the magnetic field is 
frozen into the plasma such that 



—+v-VA = Q. 

dt 


(15.6) 


This is the change of A seen by a moving plasma element. 
If, at the separatrix, this change is the same as the one 
defined by (5), plasma does not cross the separatrix, and 
magnetic reconnection does not occur. 

For reconnection, the surface and plasma velocities 
hence satisfy, on the separatrix. 


fi? -fi? - f)- 

^ ' IsiSparatriK ^ ^ ^ Ix point 
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(15T) 


(15,2) 


We can insert this result into the energy balance to yield 


where ^ 0 only inside the red -shaded region in 
Figure 15.1, implying that the first term on the right-hand 
side (RHS) of (1) can be ignored. Assuming a strictly 


^\dV 

dti 

= jdS 


p 

p — h h— ^ 

2 7-1 


v' B 


K - - 

j « ^ P 

2 2p^ l2/j„ 


(15.8) 


Figure 15.1 Schematic of the reconnection geometry. 


Equation (15.8) describes two contributions to the 
change of overall energy content of the volume enclosed 
by the separatrix: a compression or expansion contribu- 
tion (second term on RHS) and a net energy transport 
facilitated by the reconnection electric field. 

Similar arguments are readily made for mass and 
momentum transport. The analysis presented here shows 
that reconnection is a key transport mechanism, which, 
by opening otherwise closed magnetic surfaces, allows 
for significant changes to energy content. This feature is 
a critically important provider of energy sources for 
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the imiltitude of physical processes occiirriiig inside the 
reconnected plasma. 

In the following sections^ we focus on the kinetic essen- 
tials, which render reconnection possible in collisionless 
plasmas. Arguably, the most important aspect here is the 
physics of the diffusion region, which supports the electric 
field that regulates the overall transport. We will present 
aspects of the collisionless physics of the diffusion region 
for two different environments: symmetric and asym- 
metric geometries. 


of the ion Alfven velocity based on the 

reference magnitudes of magnetic field and density. The 


electric field is measured in units of and the 

pressure in units of p^- . The above length scale 


15.3. TIME-DEPENDENT RECONNECTION 


AQ5 


Time dependence and plasmoid formation have 
attracted considerable attention [e.g., Daughton et 
2006, 2009; Drake et aL^ 2006; Karirmbadl et at, 2007]. 
We here perform a detailed analysis of the correlation 
between diffusion region dimensions and reconnection 
rate. For this purpose, we employ open boundary 
condition calculations of a continuously driven recon- 
necting system. The employed code is a 2.5D variant of 
our proven 3D particle- in-cell code [Hesse et a/., 1999], 
where periodic boundary conditions have been replaced 
by open boundaries, where normal derivatives of density, 
velocity, and isotropic pressure are assumed to vanish. 
Here and in Chapter 1 6 we normalize densities by a typ- 
ical density in the current sheet, the magnetic field by 
the asymptotic value of the in-plane magnetic field. 


and Ampere's law imply that the current density is nor- 
malized to 

The initial condition consists of a poloidal magnetic 
field, a modified Harris sheet [Harris y 1962] with a current 
sheet half width of half the ion inertial length. The initial 
magnetic field is coplanar, i.e,, does not contain a guide 
field component. In a coordinate system where the x axis is 
in the direction of the initial magnetic field, the y axis is in 
the initial current direction, and the z direction completes a 
right-handed coordinate system, it is of the following form: 


Ions are assumed to be protons (mass throughout. 


and length scales by the ion ine rtial length c/ca^ where 
the ion plasma frequency ~ is evaluated 

for the reference density. Velocities are measured in units 


B and / currant denaity, t = 100.00 
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Figure 1 5.2 Two times of the simulation showing the common presence of magnetic islands at different times of 
the simulation. Note that the island shown in the bottom panel formed newly and dtd not evolve from the island 
shown in the top panel. 
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B^ = tanh(2z) + a{,^r/L^ cos(2;rx/L^ )sin(:?rz/L^ ), (15.9a) 
B^ = -afj2jitL^sm(^2jrx/L^^cos(^iiz/L^y (1 5.9b) 


The perturbation amplitude £ijj”0.1 leads to an initial 
value of the normal magnetic field of about 3% of 
The system size is L^= 100 and L^”25. The ion-electron 
mass ratio is chosen to be 25. A total of 2 x 10® 

macro particles are moved on a 1 600 x 800 grid, with an 
electron/ion temperature ratio of TJT^~ 0:2. The system 
is, at its upper boundary, subjected to a constant driving 


electric field of = OT 85. Particles are fed into the calcu- 


lation to maintain, at the upper and lower boundaries, a 
density of n^~0.2 throughout the calculation and 
removed as they leave the boundaries at 0 and x = 100. 
The calculation extends for 1 000 Alfv’en times. Figure 15.2 
displays two typical times of the simulation. 
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111 order to assess the effect of island formation on 
magnetic reconnection, we need to define the dimensions 
of the diffusion region. The electric field is given by 

£ + v x£ = — = (15.10> 

' en^ ‘ e dt ' 

Here the index e denotes electron quantities. At each X 
point, the electric field is given entirely by the non ideal 
terms on the RHS of (10). When deviating either hori- 
zontally or vertically from the diffusion region location, 
the contribution to the electric field by the convective 
term (second term on the left-hand side (LHS)) begins to 
increase, whereas the nonideal contribution decreases. 
The horizontal and vertical edges of the diffusion region 
are now defined to be the location of crossover between 
ideal and nonideal terms, i.e., 

(15.1 1> 




Equation 15.11 defines one location each to the right 
and left and top and bottom of the X point. The distance 
between the horizontal crossover points defines the length 

and the corresponding vertical distance 2^ is defined 
similarly. We perform this analysis for the dominant 
X point for each inverse ion cyclotron period during the 
simulation. The dominant X point is defined to be the 
one whose associated separatrix encloses the larges 
amount of magnetic flux. 

We expect a dependence of the reconnection electric 
field on the ratio of the two quantities. This expectation 
is demonstrated in Figure 15.3, which plots the instanta- 
neous electric field at the dominant X point against the 
diffusion region aspect ratio. 

The figure shows a clear correlation of the expected 
kind. It is noteworthy, albeit not surprising, that the 
electric field values fluctuate about the applied electric 
field. The mean reconnection electric field strength equals 
the driving value; however, reconnection remains fast at 
all times. Applying a lower or larger driving electric field 
will likely lead to a smaller or larger, respectively, average 
reconnection rate and shift the distribution accordingly. 

The point we want to make here is that the variation of 
the electric field is essentially due only to the variation of 
the diffusion region length in the x direction. For this 
purpose. Figure 15.4 illustrates the dependence of the 
electric field on the inverse diffusion region length in x, 
and Figure 15.5 depicts a histogram of encountered scale 
length in the 2 direction. 

The clear correlation of the reconnection electric field 
with the inverse diffusion region length in the x direction 
combined with the very limited variation of the diffu- 
sion region thickness, implies the variation of the 


Figure 153 Scatter plot of the instantaneous electric field at 
the major X point against the aspect ratio of the diffusion 
region. The outliers in the reconnection rates are caused by 
unusually strongly fluctuating electric fields and plasma prop- 
erties, which render an analysis very difficult. 
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Figure 15.4 Variation of the reconnection electric field at the 
dominant X point with inverse diffusion region scale length in 
the X direction. The distribution of inverse diffusion region 
lengths has a mean of 0.54 and a standard deviation of 0.27. 
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Figure 15.5 Histogram of diffusion region thicknesses 
throughout the calculation. The figure demonstrates only very 
small variations of this parameter. The distribution has a mean 
of 0.67 with a standard deviation of 0.09. 


reconnection electric field is primarily due to variations 
of the diffusion region dimensions in the x direction. 
This latter variation is caused by changes to the electron 
magnetization in the horizontal direction, which extends 
or contracts the extent of the diffusion region. Electron 
magnetization is primarily affected by the strength of 
the magnetic field. 

It appears evident that the growth of tearing mode 
islands, which form in the reconnecting ctirrent layer, are 
the prime cause of the variation of this magnetization. As a 
tearing mode island becomes nonlinear, it will begin to 
magnetize electrons and thus shorten the electron bounce 
width in the X direction. This leads to an effective short- 
ening of which, as we have seen above, will typically lead 
to a faster reconnection rate. Therefore, we conclude that 
island formation in the current layer associated with an 
active X point can increase the rate of magnetic recon- 
nection. However, we note that reconnection can be fast 
for longer diffusion regions, presumably related to the 
absence of islands in the proximity of the X point under 
consideration. Figure 1 5.4 clearly shows that the rate virtu- 
ally never drops below £=0.12. We point out that the 
applied electric field exceeds the empirical minimum values. 
Hence it is conceivable that plasmoid formation may be 
necessaiy to raise the average reconnection rate to that 
required for the flux transport prescribed by the external 
driving. 


15 . 4 . RECONNECTION LINE ORIENTATION FOR 
ASYMMETRIC RECONNECTION 

In this section, we review recent results [Hesse et aL, 

2013] pertaining to another aspect of magnetic reconnec- 
tion: merging of plasmas from asymmetric inflow regions. 

Such asymmetry is commonly found at the magneto- 
pause of Earth as well as on other planets. Typically, the 
shear angle of the magnetic field is different from the 1 80® 
modeled above, which, in comparison, is represented by 
the presence of a guide magnetic field. Here a key 
question is which direction a reconnection line will follow: 

Will it be normal to the plane wherein the guide field is 
constant [Sonnerup, 1974] or will it choose a different 
direction [Gonzales and Mozer, 1974; Cowk}\ 1976]? In 
order to address this question, we model a sample system 
with parameters set up similar to Pritchett [2008]. The 
polo id al magnetic field is of the form 

B^ = i^5^tanh[zll\ (15J2) 

with a guide magnetic field 

= (15.13) 

The initial ion and electron densities are chosen to 
exhibit asymptotic values of n-1 for z>0 and n-03SS 
for z<0\ 

n. ~ 1 /3.- [z//]/3. (15.14) 

The initial current layer half-width is /= 0.5. A temper- 

ature of ryr^=0.2 and a temperature of T= T^= 1.5 
are employed. A small, X-type, initial perturbation is AQ5 
introduced into current density and magnetic fields with 
an amplitude of 5£ = 0.1, leading, e.g., to a form of B^ 
similar to (9b). 

Particles, initialized as drifting Max well ians, are 

integrated on a grid of 1000x800 cells representing a 
system of size L^ = 64 and L^=25.6. Periodic boundary 
conditions are employed in the x direction, whereas the 
particles are specularly reflected at the upper and lower 
boundaries The ion-electron mass ratio is chosen to be 
w/w^=25. A total of 1.6x10^ macro particles are 
employed for the simulations. 

In addition to investigating the reference case above, we 
further study the guide field model in coordinate systems 
of different orientation, derived from the one above by 
rotation around the z axis. The rotation angle a is defined 
in Figure 15.6. For each simulation in a rotated coordinate 
system, the entire set of vector quantities defined above is 
transformed into the rotated coordinate system, and the 
numerical model is applied to the new initial condition. 

The overall evolution of the reference system, with uni- 
form guide field of unity (i.e., a = 0), is illustrated for 
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Figure 15.6 Graphics of the asymptotic magnetic field values 
in a coordinate system rotated by an angie a about the t axis. 
The indices u (up) and d (down) refer to asymptotic values 
above and below the current layer, respectively. The globally 
constant guide field 1 for«=0 becomes spatially dependent 
after rotation. The asymptotic values are shown on the / axis, 

three different times in Figure 15.7. The figure shows a 
transition from initially symmetric behavior to the 
development of a significant left-right asymmeti^ toward 
the end of the simulation period = 80), 

The temporal evolution of the reference reconnection 
rate is shown in Figure 1 5.8. The figure demonstrates that 
Sonnerup^s [1974] prediction of fast reconnection in the 
presence of a guide field is correct: The peak reconnec- 
tion rate is still about 0,05. The simulation never attains a 
steady state, a consequence of both periodic boundary 
conditions in the x direction and the closed top and bot- 
tom boundaries, which imply limited amount of magnetic 
flux and energy to tap into. Open boundary conditions 
for asymmetric systems are difficult to implement, which 
led us to choose the present implementation. 

A reasonable assumption has the X line naturally orient- 
ing itself in such a way to obtain the fastest reconnection 
rate. The most important question is therefore whether 
faster reconnection can be obtained if a* 0. Conducting a 
set of simulations, we find substantial variations of the 
reconnection rates depending on the orientation of the 
coordinate system. Figure 15.9 displays the time evolution 
of the reconnection rates for all runs obtained from the 
guide field reference model by coordinate system rotation 
and illustrates this strong dependence. 

Figure 1 5.9 exhibits a maximum reconnection rate for 
a = -14.87'^. Adjacent values of the rotation angle yield 
almost identical reconnection rates, indicating that the 
maximum is not strongly localized. Extending rotation 



Figure 15.7 Magnetic field (white lines) and current density 
(color) evolution for the reference run, for which the initial 
guide field is uniform and of unit value. 



Figure 15.8 Time evolution of the reconnection electric field 
for the reference run for which the initial guide field is uniform 
and of unit value. 
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Figure 15-9 Time evolution of the reconnection electric field 
for the entire set of runs derived from rotating the frame of the 
guide field calculation by an angle «. The different colors 
denote different runs, and the angles are denoted in the figure. 


Figure 15-10 Plot of peak electric field values and predictions 
based on the magnetic energy available for magnetic recon- 
nection. The magnetic-energy- based prediction exhibits an 
excellent match. 


angles more significantly to either side, however, yields 
rapidly decreasing values of the reconiiection electric 
fields. It is noteworthy that the reconnection evolution for 
a = -26.57*^, for which the reconnection magnetic fields 
are equal on both sides, still leads to considerably slower 
reconnection rates than the runs for smaller rotation 
angles. Therefore, symmetry of the in-plane magnetic 
field does not lead to the fastest reconnection rates. 

Furthermore, the reconnection rate has to vanish for 
rotations, for which one of the two rotated magnetic field 
components 

cos(a) + B^ sin(a), (15.15) 

B^, - B^ cos(a)-i- sin(a) (15.16) 

vanishesv For the system investigated here, ^^=1.5 and 
B= —0,5, We find ^ ^ = 0 for the rotation value a~ — 57^. 

d u 

In order to develop a scaling relation for the reconnec- 
tion electric field, we consider the maximum value of the 
entire time evolution for a fixed rotation angle. The key 
question is how the maximum values shown here relate to 
physical parameters. We will approach this issue from two 
different anglea First, we will make a reasonable assump- 
tion: reconnection rates should scale with the available 
magnetic energy. Because of the constancy of the total 
volume, this assumption implies that reconnection rates 


should be proportional to the magnetic energy density. 
Since the energy densities available for reconnection can 
differ across the current layer, our assumed proportion- 
ality implies 

(15J7) 

where the factor should depend on parameters 
such as density and total magnetic field, which are held 
constant in the present investigation. Figure 15.10 
demonstrates that the model (17) does an excellent job 
representing the maxima of the modeled reconnection 
rates. 

We can use our model (17) to derive the rotation angle, 
for which this energy product assumes its maximum. 
Using the rotated magnetic field expressions (15) and 
(1 6), it is a straightforward exercise to derive the equation 

B B — B ^ 

cot"a^-2 \ cota„-l = 0 (15.18) 

for the angle for which (17) attains its maximum. 
The relevant solution is 

=- 14 . 87 °, 

matching our empirical determination. 
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It is further interesting to note that equation 15, IS also 
describes exactly the angle which bisects the angle bet- 
ween the asymptotic magnetic field on both sides of the 
current layer, a result readily derived from geometric con- 
siderations, This result applies to more than the system 
modeled here and is hence quite generic: The ha If- angle 
direction maximizes the magnetic energy available for 
reconnection. 

In summary, our results indicate that the magnetic 
reconnection line in asymmetric systems is preferentially 
oriented in such a way that it bisects the direction of the 
asymptotic magnetic field direction on both infiow sides. 
This orientation is identical to the one for which the 
product of available magnetic energy is maximized. 

15 . 5 . SUMMARY 

In this chapter, we presented results pertaining to three 
different aspects of magnetic reconnection. The first 
aspect, transport of magnetic flux, mass, momentum, 
and energy, is generic to all geometries and environments 
in which reconnection is operating. Focusing on energy 
transport, we could show that magnetic reconnection 
enables effective energy transport beyond what is possible 
by compression or expansion of a volume enclosed by a 
separatrix. We showed that the critical quantity providing 
the amplitude of the transport rate is exactly the recon- 
nection electric field, measured at the X point. While not 
developed here, it is straightforward to develop similar 
analyses for mass and momentum transport. 

Having demonstrated that effective energy transport is 
impossible without magnetic reconnection, we embarked 
on a discussion of the second topic: time dependence. 
Using our open boundary condition models, we executed 
an extended calculation of a driven, continuously recon- 
necting, symmetric current sheet. 

For the purpose of the analysis of the model run, we 
defined diffusion region dimension through the region, 
where nonideal contributions to the electric field exceed 
the electric field contribution by electron convection. The 
magnitude crossover points were identified with the 
horizontal and vertical diffusion region edges. We devel- 
oped an automated algorithm to calculate, for each 
output time step, these dimensions and take note of the 
instantaneous value of the reconnection electric field. 

We then used the results for a statistical correlation. 
Not surprisingly, we found that rates tend to correlate 
well with the aspect ratio of the diffusion region. More 
surprisingly, we discovered that the aspect ratio variation 
was almost exclusively due to variation of diffusion 
region lengths, whereas diffusion region width showed 
relatively limited variationa We explained this phenomenon 
by the repetitive formation and grov^^h of plasmoids, which, 
when grown past the linear stage, begin to magnetize 


electrons. This magnetization changes the crossover loca- 
tion between nonideal and convection electric fields and 
renders the diffusion region shorter. Thus, the aspect 
ratio becomes larger, and, on average, the reconnection 
electric field grows. This modulation is of the order of 
50%, implying that reconnection is fast even in the 
temporary absence of magnetic islands. Therefore, our 
conclusion is that island formation is not necessary for 
fast reconnection, but it typically serves to enhance 
the reconnection rate substantially and in a very time- 
dependent manner. 

As the third part of this chapter, we briefly reviewed 
recent results pertaining to the direction of the reconnec- 
tion line in aS 5 rmmetric magnetic reconnection. Starting 
with an asymmetric configuration with a guide magnetic 
field, which, in many ways, is typical for reconnection 
geometries expected at a planetary magnetopause, we 
performed a set of kinetic simulations of this system in 
rotated coordinate systems. We found a surprising varia- 
tion in reconnection rates, even more surprising since nei- 
ther constant guide field nor symmetric in*plane magnetic 
fields led to the fastest reconnection rate. Instead, we 
found empirically that the reconnection rates scaled well 
with the product of energy densities of both inflow 
regions. The direction which maximizes this product 
turned out to be also halving the angle between the 
asymptotic magnetic field directions on both sides of the 
current layer. Clearly, this result requires further testing, 
including for different density asymmetries but it is highly 
suggestive that the X line will orient itself such that it 
bisects the directions on the magnetic fields in the two 
inflow regions. 

Finally, one might ask the question of how these results 
pertain to the magnetospheres of planets other than 
Earth. The general results presented in section 15.2. are 
universal, i.e,, they apply, at least in principle^ also to 
planetary magnetospheres. As the analysis in subsequent 
sections has been conducted in dimensionless units, a 
large variety of applications can be developed by employ- 
ing appropriate dimensions in the interpretation of the 
results. This freedom is an application of the 'Wlasov 
scaling” (Schindler, 1969], whereby Vlasov scaling laws 
permit rescaling of solution to different applications. 
Such rescaling will greatly extend the applicability of the 
results presented here; however, just as in the terrestrial 
magnetosphere, rescaling will fall short of covering all 
relevant scenarios. Therefore, much more research needs 
to be conducted, and many more exciting and surprising 
discoveries lie in wait. 
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